Strange attractor for the renormalization flow for invariant tori of Hamiltonian 

systems with two generic frequencies 

C. Chandre and H.R. Jauslin 
Laboratoire de Physique-CNRS, Universite de Bourgogne, B.P. 47 870, F-21078 Dijon, France 



0\' 



m 



>■ 

<n: 

o: 
o. 
\o- 

o: 

a^ 
a^ 

I 

o 



% 



We analyze the stability of invariant tori for Hamiltonian 
systems with two degrees of freedom by constructing a trans- 
formation that combines Kolmogorov-Arnold-Moser theory 
and renormalization-group techniques. This transformation 
is based on the continued fraction expansion of the frequency 
of the torus. We apply this transformation numerically for 
arbitrary frequencies that contain bounded entries in the con- 
tinued fraction expansion. We give a global picture of renor- 
malization flow for the stability of invariant tori, and we show 
that the properties of critical (and near critical) tori can be 
obtained by analyzing renormalization dynamics around a sin- 
gle hyperbolic strange attractor. We compute the fractal dia- 
gram, i.e. the critical coupling as a function of the frequencies, 
associated with a given one-parameter family. 

PACS numbers: 05.45.Ac, 05.10.Cc, 45.20.Jj 



I. INTRODUCTION 

Invariant tori play a fundamental role in the stabil- 
ity of Hamiltonian systems. For two degrees of freedom, 
they act as barriers in phase space. Renormalization- 
group (RG) transformations have been proposed to un- 
derstand the breakup of invariant tori for area-preserving 
maps y,0] and for Hamiltonian systems with two degrees 
of freedom S-S]. The RG-transformations described in 
Refs. |§,0,|| combine an elimination of the irrelevant 
part of the perturbation, and a rescaling of phase space 
which is adapted to the frequency of the specific invariant 
torus. These transformations are constructed as canoni- 
cal changes of coordinates. They have been implemented 
numerically for the golden mean torus (with frequency 
7 = (V5 — l)/2). Based on the numerical analysis of 
the renormalization flow, the following picture emerges : 
A trivial flxed point Hq (which is an integrable Hamil- 
tonian) of the RG-transformation characterizes the do- 
main where the Hamiltonians have a smooth invariant 
torus with frequency 7, i.e., all Hamiltonians attracted 
to Hq by renormalization are locally canonically conju- 
gated to an integrable Hamiltonian. The boundary of 
the domain of attraction of Hq is the domain where the 
tori are critical, at the threshold of their breakup. This 
boundary (named the critical surface) is expected to be 
of codimension 1, and to coincide with the set of Hamil- 
tonians with critical coupling, beyond which they do no 
longer have an invariant torus with frequency 7. The 
critical surface is expected to be the stable manifold of a 
nontrivial fixed point H^. This global picture is still at 



the stage of conjecture : In the perturbative regime, it 
is supported by rigorous results Q . From the numerical 
analysis of the neighborhood of H^, , the scaling properties 
of critical tori PJTQ] are found by these renormalization 
transformations : It is a strong argument in favor of the 
link between H^ and critical invariant tori. Furthermore, 
this renormalization approach gives a method to calcu- 
late the critical coupling for a given one-parameter family 
(the value of the parameter for which the considered in- 
variant torus is critical). It has been verified that this 
critical coupling coincides with the critical couplings ob- 
tained by other methods like Greene's criterion |11[| or 
Laskar's frequency analysis p^ . 

Since the relevant quantity for stability is the ratio 
between the two components of the frequency vector, 
we consider an invariant torus with frequency vector 
a;o = (w,— 1) where the frequency cu is irrational (the 
frequency vector must satisfy a Diophantine condition). 
The purpose of the RG-transformations we define is to 
analyze the stability of this invariant torus for Hamilto- 
nians with two degrees of freedom. The idea is to set 
up a transformation 7?, as a canonical change of coordi- 
nates that maps a Hamiltonian H into a rescaled one 
TZ{H), such that irrelevant degrees of freedom are elimi- 
nated. The transformation TZ should have the following 
properties : TZ has an attractive fixed set (trivial fixed 
set) of integrable Hamiltonians. Every Hamiltonian in 
its domain of attraction V has a smooth invariant torus 
with frequency vector Wq. The aim is to show that there 
is another fixed set which lies on the boundary dV (the 
critical surface) and that is attractive for every Hamilto- 
nian on dV. 

Renormalization-group approaches for the breakup of in- 
variant tori for arbitrary frequency have been proposed 
for circle maps |13-1^] and for area-preserving maps [[l8| . 
The boundary of Siegel disks has also been analyzed 
by renormalization |19| . The numerical results suggest 
that statistical self-similarity characterizes the considered 
problems at criticality. It has been conjectured that it 
can be described by an ergodic attractor of a renormal- 
ization transformation. 

The numerical implementation of the RG-transformation 
for Hamiltonian systems with two degrees of freedom, 
gives support to this picture : We conjecture that the 
properties (scaling factors and critical exponents) of crit- 
ical tori can be obtained by analyzing renormalization dy- 
namics around a single strange chaotic attractor. Each 
critical torus displays a sequence of critical exponents, 
which are the eigenvalues of the linearized renormaliza- 



tion map along its trajectories. Two different critical 
tori display the same set of critical exponents but with a 
different probability distribution. In that sense, we can 
speak of a single universality class describing critical in- 
variant tori for Hamiltonian systems with two degrees of 
freedom. 

In Sec. y, we construct the RG-transformations for 
an invariant torus with bounded entries in the continued 
fraction expansion of its frequency. 



Ill 



we ap- 



In Sec. 

ply this construction numerically for a set of frequencies 
which have only 1 and 2 in their continued fraction expan- 
sion. We show that a critical strange attractor can be ex- 
pected to describe the properties of critical invariant tori. 
In Sec. [V, we compute a fractal diagram which is the 
set of critical couplings Sc{uj) for a given one-parameter 
family of Hamiltonians. In Sec. [V|, we construct a simple 
approximate RG-transformation that gives qualitatively 
all the relevant features of the renormalization dynamics. 



the torus with frequency vector ljq — (w, — 1). The next 
smaller scale is characterized by uf^^ = (l/ag,— 1). The 
renormalization transformation, denoted TZao-, changes 
the coordinates such that the next smaller scale becomes 
the main one, i.e., the main scale in the new coordinates 
is characterized by uj'^ = oj^,-^ . It acts as a microscope 
in phase space : it looks the system at a smaller scale in 
phase space, and at a longer time scale. 

We consider Hamiltonians H written in terms of ac- 



tions A = {Ai. 
of the form 



A2) e R^ and angles V = (vJi, V'2) G T^ 



H{A,^) = H^{A) + V{n ■ A,^), 



(2.2) 



where 17 = (1, a) is a vector not parallel to the frequency 
vector Wo, and Hq is given by 

Ha{A)^u,o-A+\{n-A)\ 



II. RENORMALIZATION TRANSFORMATION 

We describe the renormalization scheme for a torus 
with arbitrary frequency vector Wq = (w,— 1) where 
uj G]0, 1[. This renormalization relics upon the contin- 
ued fraction expansion of ld : 



1 



[ao,ai,. 



ao 



ai 



The best rational approximants of w are given by the 
truncations of this expansion : pk/lk — [ooi ^ii ■ ■ • 1 Ofe = 
00]. The corresponding periodic orbits with frequency 
vectors ijJ,^^. = (j>k/Qk,—i) accumulate at the invariant 
torus. We call "resonance", a vector v>k = {qk,Pk) (for 
fc > 1) which is orthogonal to w^^,. The word resonance 
refers to the fact that the small denominators ujq ■ v^ that 
appear in the perturbation expansion are the smallest 
ones, i.e., |a;o • i>'A;| < l^^o • i'\ for any v = {q,p) different 
from zero and i/^, and such that \q\ < qk+i PQl. The 
sequence of resonances satisfies Jwq • I'k+i \ < \<^o-i^k\ and 
linifc^oo |'*^o ■ i^fcl = Oj as can be seen from the relation 



i^k = Nao ■■■Na^_^ua, 
where Uq — (1, 0) and Na^ denotes the matrix 



(2.1) 



Na, 



ai 1 
1 



Moreover, uiq • f fc and Wq • f fc+i are of opposite sign (as 
the stable eigenvalue of N^i is negative) ; thus the torus is 
approached from above and from below by the sequence 
of periodic orbits with frequency vectors {a;^,^}. 

The main scale of the torus is characterized by the 
periodic orbit with frequency vector w^^ = (0,-1) for 



We notice that the invariant torus with frequency vector 
Wo is located, for i?o, at A such that Vt, ■ A — and 
ljq ■ A = E, where E is the total energy of the system. 
The RG-transformation TZao consists of four steps : 
(1) A shift of the resonances constructed from the con- 
dition f 1 I— » f : we require that cos[(ao, V) ■ ^p] — 
cos[(l, 0) • if']. This change is done via a linear canonical 
transformation 

{A,^)^{A\^')^{N-^lA,Na,v). 

This step changes the frequency vector ujq into ui'q = 
(cj',— 1), since NagU^o = —uju'q, where the frequency uj' 
is given by the Gauss map 



1-^ Lj' = UJ — [tJ ] 



(2.3) 



and [i^^^] denotes the integer part of uj^^. Expressed 
in terms of the continued fraction expansion of the fre- 
quency, it corresponds to a shift to the left 



UJ = [ao,ai,a2. 



uj' = [ai,a2,a3, . 



The sequence of the resonances {i^k} is mapped into the 
sequence 

The linear term in the actions of Hq is changed into 
—ujlJq ■ A'. 

(2) We rescale the energy by a factor oj^^ (or equivalently 
time by a factor lu), and we change the sign of both phase 
space coordinates {A,ip) 1-^ {—A,—(p), in order to have 
uJq as the new frequency vector, i.e., the linear term in 
the actions of Hq is Wg • A. Furthermore, il = (l,a) is 
changed into J7' = {l,a') = (1, (oq + a)^^)- The map 
a I— > (ao + a)^^ is the inverse of the Gauss map, i.e. if 
a — [60 , ^1 , . . .] , then a' — [aQ,bo,bi, . . .]. We remark that 



if a has the continued fraction expansion a — [feg, 61, • • ■], 
and if we define the two-sided sequence 

[a|cj] = [. . .,62,6i,feo|ao,ai,a2.. .], 

then the map \o\lo\ i— > [a'Iti''] corresponds to the (two- 
sided) Bernoulh shift 



, 62,61, 6o|ao,ai,a2, . 



■ ,62,6i,6o,ao|ai,a2. 



(3) Then we perform a rescahng of the actions : H is 
changed into 



A 



H'{A,^)^\H[-,^ 



with A = \{H) such that the mean value (i.e., the average 
over the angles) of the quadratic term in the actions in H' 
is equal to 1/2. This normalization condition is essential 
for the convergence of the transformation. After Steps 1, 
2 and 3, the Hamiltonian expressed in the new variables 
is 



H'{A,ip) = \uj-^H [ --NA,-N-^^ 



For H given by Eq. (|2.2D, this expression becomes 



(2.4) 



H'{A, cp) = a;[, • A + ^^(ao + af{^' ■ Af 



+\uj-^V 



0-13 + a , 



X 



n' ■ A,-N-^ip] . (2.5) 



Thus the choice of the rescaling in the actions (Step 3) is 



X = Lu-\ao + a)^{l + 2{V''^^)), 



(2.6) 



where {V^^') denotes the mean value of the quadratic 
part of V, in the (Jl' • A)-variable. 

(4) The last step is a canonical transformation that elim- 
inates the nonresonant part of the perturbation in H'. 
The choice of which part of the perturbation is to be con- 
sidered resonant or not is somewhat arbitrary. The set 
of nonresonant modes includes the modes of the pertur- 
bation which are sufficiently far from the resonances in 
order to avoid small denominator problems during this 
elimination step. A reasonable choice for the nonreso- 
nant modes is the set /~ of integer vectors 1/ e Z^ such 
that 1 1^2 1 > 1 1^1 1- A mode which is not an ele ment of /^, 
will be called resonant. We notice that Eq. (2.1) defin- 
ing i/k = {qk,Pk) shows that qk > Pk for fc > 0, i.e., the 
resonances are not elements of /^ . From the form of the 
eigenvectors of Na^ , one can see that every v ^1?\ {0} 
goes into /" after sufhciently many iterations of matri- 
ces Nai (as the eigenvector of N~^ with an eigenvalue 
of norm larger than 1 points into /^). In other terms, a 
resonant mode at some scale turns out to be nonresonant 
at a sufficiently smaller scale. We notice that is not an 
element of /^, i.e., it is resonant. 



We eliminate completely all the nonresonant modes of the 
perturbation by a canonical transformation, connected 
to the identity, which is defined by iterating KAM-type 
transformations (one iteration reduces the nonresonant 
modes of the perturbation from e to e^). 
One iteration step is performed by a Lie transformation 
Us ■ {A, if) 1-^ {A\ if') generated by a function S{A, ip). 
The expression of the Hamiltonian in the new coordinates 
is given by 



H' = HoUs = e^H 

= H + {S,H] + ^^{S,{S,H}} 



(2.7) 



where { , } is the Poisson bracket of two functions of the 
action and angle coordinates 

Sf \^ ^ ^_^ ^ 
^•''^^ dif'dA dA' dip' 

and the operator 5* is defined as SH — {S,H}. The 
generating function S is chosen such that the order e of 
the nonresonant part (the modes in /^) of the perturba- 
tion vanishes. We construct recursively the Hamiltonians 
Hn, starting with Hi = H', such that the limit iJoo is 
canonically conjugated with H' but does not contain non- 
resonant modes. One step of this elimination procedure, 
Hn I— > Hn+i, is done by applying a change of coordi- 
nates Un such that the order of the nonresonant modes 
of Hn+i = Hn o Un IS £^, whcrc e„ denotes the order 
of the nonresonant modes of i/„. At the n-th step, the 



order of the nonresonant modes of Hn is s, 



2"- 



where Sq 



is the order of the nonresonant modes of H' . If this pro- 
cedure converges, it defines a canonical transformation 
Us = Ui 0U2 ° ■ ■ ■ ° Un o • • • , such that the final Hamil- 
tonian Hrx, ^ H' o Us does not contain any nonresonant 
mode. 

The specific implementation of this step can be per- 
formed in two versions : the first one is a transformation 
acting in a space of quadratic Hamiltonians in the actions 
of the form 

H{A, <p)=uJo-A+ ^m{v){n ■ A)2 + g(^)n . A + f{<p), 

(2.8) 

following Thirring's approach [£l|j^,|2^, where m, g, and 
/ are three scalar functions of the angles. The second ver- 
sion is a transformation acting on a more general family 
of Hamiltonians, with a power series expansion in the 
actions. 



HiA, <p)^u:o-A + J2 f^'^ i<P)i^ ■ ^)^ 



(2.9) 



i=o 



following Ref. B, where {/'-'■'} are scalar functions 
of the angles. The key-point for the implementa- 
tion for quadratic Hamiltonians is that the KAM- 
transformations do not reduce the nonresonant modes of 



the quadratic part of the Hamiltonian. These transfor- 
mations are performed by Lie transformations generated 
by functions linear in the actions. Following this proce- 
dure, the iterated Hamiltonians iJ„ remain quadratic in 
the actions. 

The eq uat ions defining the elimination part for Hamilto- 
nians (2.8), are given in Refs. ||7|,p3|]. Concerning those 
for Hamiltonians ( |2.9| ), they are given in the Appendix. 

The convergence of the elimination procedure (step 4, 
definition oiUs) has been rigorously analyzed in Ref. 0] 
for the second version of the transformation. It has been 
proven that for a sufficiently small perturbation, the elim- 
ination converges. Concerning the quadratic case, we 
lack, at this moment, a theoretical framework to prove an 
analogous theorem. The convergence of the elimination 
procedure outside the perturbative regime is observed in 
both cases numerically. 

In summary, the RG-transformations for a given torus, 
act as follows : First, some of the resonant modes of the 
perturbation are turned into nonresonant modes by a 
frequency shift and a rescaling in phase space. Then, 
a KAM-type iteration eliminates these non-resonant 
modes, while producing some new resonant modes. 
Essentially, both versions of the RG-transformation give 
qualitatively and quantitatively the same results for the 
class of invariant tori considered in this article. The 
main advantage of the first version is that it is numeri- 
cally more efficient as one only works with three scalar 
functions of the angles. But from a theoretical point of 
view, there are some advantages to work with the second 
version, since the first one leads asymptotically to non- 
analytic Hamiltonians for the attracting fixed sets (in the 
quadratic part in the actions) (t). 



III. NUMERICAL RESULTS 

For a given frequency u (with bounded entries in its 
continued fraction expansion), the numerical implemen- 
tation of the RG-transformation shows that there are 
two main domains in the space of Hamiltonians : one 
where the iteration converges towards a family of in- 
tegrable Hamiltonians (trivial fixed set), and the other 
where it diverges to infinity. These domains are sepa- 
rated by a surface called critical surface (which is con- 
jectured to coincide with the set of Hamiltonians which 
have a non-smooth critical invariant torus of the consid- 
ered frequency) and denoted S{uj) in what follows. 
The domain of attraction of the trivial fixed set is the 
domain where the perturbation of the iterated Hamilto- 
nians tends to zero. However, the renormalization trajec- 
tories do not, in general, converge to a fixed Hamiltonian 
but to a trajectory related to the Gauss map (2^). The 
trivial fixed set is composed of Hamiltonians of the form 



where w; = (w;,— 1) and fli = (1,0;;). The renormaliza- 
tion map transforms the vectors uji and il; following the 
Gauss map 

ai+i = {ai + [^r^])~^- 

Thus the trivial fixed set can be a fixed point, a periodic 
cycle, or in general a set of Hamiltonians labeled by a 
trajectory of the Gauss map, i.e., by the asymptotic en- 
tries of the continued fraction expansion of to. 
Correspondingly, on the critical surface 5(a;), the renor- 
malization flow converges to a periodic cycle if the fre- 
quency a; is a quadratic irrational (its continued fraction 
expansion is asymptotically periodic), and to a low di- 
mensional attracting set which is a strange chaotic at- 
tractor, in the case where uj is not quadratic. The non- 
trivial attractor has a codimension 1 stable manifold, i.e., 
one expansive direction transverse to the critical surface 
S{lu). 

We take the definitions as formulated by Grebogi et 
al. p3] : An attractor of a map is called strange if it 
is not a finite number of points, nor a piecewise differen- 
tiable set. An attractor is chaotic if typical orbits on it 
have a positive Lyapunov exponent. 

Quadratic irrational frequencies- We start by analyz- 
ing the effect on lo and a, of s renormalization steps. 
Denote by {5^} the continued fraction expansion of a : 
a= [60,^1, •• ■]■ The renormalization 7?.a^_j7?.a3 _ 2 ' ' '"^ao 



into [a 



S,"5 + l, ■ 



.], and a into 



Hi{A)=LJi-A+]-{^i-Af, 



(3.1) 



changes uj — [ao,ai, 

[as_i,as_2,.. .,ao,&o,fei, • • ■]■ 

If UJ has a periodic continued fraction expansion of pe- 
riod s > 1, i.e., oj = [(ai, . . . , as)oo], we notice that a 
converges to [(a^, . . . ,ai)oo]- Therefore Q, converges to 
the unstable eigenvector of the matrix Na^ ■ ■ ■ Na^ (the 
stable eigenvector of this matrix is Wq)- 
In that case, the RG-transformation has no fixed points 
but two periodic cycles with period s : a trivial cycle 
which is attractive in all the directions in the space of 
Hamiltonians, and a hyperbolic (critical) cycle with a 
codimension 1 stable manifold. The trivial cycle charac- 
terizes the domain of Hamiltonians that have a smooth 
invariant torus of the considered frequency, and the non- 
trivial cycle, the domain where the torus is critical. 
These periodic cycles can be equivalently considered as 
fixed points of the RG-transformation defined by the 
composed operator TZg = TZa^'R-a^_i ■ ■ ■ Tlai ■ 
The nontrivial fixed point of TZg associated with u de- 
fines a universality class that we characterize with criti- 
cal exponents such as the total rescaling of phase space 
(product of the s rescalings) , and the unstable eigenvalue 
of the linearized map around the fixed point. 
The interpretation of these attracting cycles (trivial and 
non-trivial) of the RG-trajectories in terms of the struc- 
ture in phase space, of the invariant tori is the following : 
In the perturbative regime, there exists a geometrical ac- 
cumulation of a sequence of periodic orbits. The fact that 



it also happens in the critical case, but with a nontrivial 
ratio, implies universal self-similar properties of the crit- 
ical torus [10. 

Two frequencies uji and cj2 having the same periodic tail 
(and different first entries) in their continued fraction ex- 
pansions (wi and LU2 are called equivalent), belong to the 
same universality class : After a finite number of steps of 
the RG-transformation, the equations defining the renor- 
malization trajectories for cui are the same as those for 
LU2. Thus they have the same critical exponents and scal- 
ing factors. The initial integers in the continued fraction 
expansion are irrelevant for the breakup of the torus. 
Associated with this nontrivial fixed point, we also have 
nontrivial fixed sets related to the nontrivial fixed point 
by symmetries ^,^,H,g3[7 which therefore belong to the 
same universality class. 

Nonquadratic irra- 

tional frequencies- The RG-transformations constructed 
in Sec. O are in principle defined for an arbitrary irra- 
tional frequency oj €]0, 1[. They are based on the con- 
tinued fraction expansion of the frequency. However, the 
numerical analysis of the transformation on the critical 
surface 5(w), requires to have a large number of entries 
in the continued fraction of uj (in order to iterate the 
RG-transformation on Siuj) a sufficient number of times 
so as to observe the convergence to an attracting set). 
We consider frequencies whose entries in the continued 
fraction expansion are chosen randomly in a finite set of 
integers, for instance, 1 or 2. Each entry in the continued 
fraction expansion is chosen independently of the others 
according to the probability P(l) = p and P(2) = 1 — p. 
For p = 1, the frequency is equal to the golden mean 7, 
and it is equal to \/2 — 1 for p = 0. 

For a given frequency uj with p G [0,1], the numerical 
analysis shows that there are two main attracting sets in 
the space of Hamiltonians : a trivial fixed set composed 
by integrable Hamiltonians, and a nontrivial one which 
lies on the critical surface 5(w). 

Figure |l| depicts a projection of the critical attractor ob- 
tained for a typical frequency with p — 1/2, and Fig. g 

for a frequency with p = 9/10, on the plane (/ii , /ii ) 
where f^i denotes the Fourier coefficient Ui of the func- 



tion /('^ of a Hamiltonian ( |2.9| ) on the critical attrac- 
tor. For numerical purposes, these pictures have been ob- 
tained by truncating the Hamiltonians on which the RG- 
transformation acts, in power series (we neglect higher 
orders than (Jl • A)^^ , with J = 5), and in Fourier series 
of the functions f^^^^ (wc neglect the Fourier modes with 
maxj ji^ij > L, with L = 5). 

A RG-trajectory on the critical attractor displays dif- 
ferent scaling factors and critical exponents at each point 
of the trajectory. The distribution of these exponents de- 
pends only on the frequency of the torus. 
In Fig. H, we depict the values of the rescahngs (A,.. Ai+i ), 
where Xi denotes the value of the rescaling value (2.6) af- 



ter i iterations on the critical attractor. This picture is 
obtained for a frequency lu with p — 1/2, and the cut-off 
parameters are J — 5 and L = 5. There are four main 
parts in this graph. They correspond to the four possible 
changes of the first entry (to the second one) of the con- 
tinued fraction of the frequency : 1^1, 1^2, 2^1, 
2^2. There are two values for which A^+i = Xi : 4.339 
and 14.871 which are the values for the rescaling for the 
frequencies (\/5 - l)/2 and \/2 - 1 |lO|J2|]. 



If we consider two frequencies associated with a same 
probability P(l) — p, the numerically computed projec- 
tions on the Fourier modes of the perturbation, of the 
critical attractors for each frequency, look identical. This 
suggests that the properties of critical tori do not depend 
on the order of the entries in the continued fraction ex- 
pansion, but only on the distribution of the values of the 
entries. In other words, critical tori with frequencies uji 
and 0^2 associated with a same probability distribution 
in the continued fraction expansion, are statistically self- 
similar. We notice that this discussion encompasses also 
the case of two equivalent frequencies. 

For two frequencies associated to different distributions 
of the entries in the continued fraction expansion (for 
different values of p), the support of projections (on the 
Fourier modes of the perturbation) of the critical attrac- 
tors are the same, but the renormalization trajectories 
visit the attractor with different densities. 
For instance, concerning the projection of the attractor 
found for a frequency with probability p = 9/10 in Fig. g, 
since p > 1/2 the renormalization trajectory is more con- 
centrated around the fixed point (or around a periodic 
cycle related to this fixed point by symmetry) obtained 
for the golden mean (p = 1). 

In order to describe a more precise image of the renor- 
malization flow, we consider an enlarged space where 
we add to the space of Hamiltonians a direction corre- 
sponding to the frequency of the considered torus. The 
RG-transformations act in this enlarged space consisting 
of couples {H, to) (we iterate a renormalization operator 
for a given Hamiltonian H and for a given torus charac- 
terized by its frequency lu). The numerical results lead 
to the following conjecture : The enlarged space is di- 
vided into two main parts : one where the RG-iteration 
converges to a trivial fixed set, that consists of couples 
{Hi,uji), where Hi is integrable (of the form (3.1)). This 
trivial fixed set attracts all points {H, uj) such that H has 
a smooth invariant torus with frequency lu. The dynam- 
ics on it is determined by the Gauss map. The domain of 
attraction of this trivial fixed set is bounded by a critical 
surface S which is the set of the critical surfaces S{uj), 
i.e., composed by points {H,u;) such that H has a non- 
smooth invariant torus with frequency lu. The surface S 
has a fractal structure : In Sec. |^, we display a section 
of this surface, and analyze the fractal diagram of a given 
one-parameter family. On the critical surface iS, there is 



a single critical attractor to which all Hamiltonians on S 
are attracted. On these attractors (trivial or critical), we 
have subsets that are periodic cycles of all the periods; 
these cycles correspond to quadratic irrational frequen- 
cies. For instance, we have two fixed points on the non- 
trivial attractor : one corresponding to the golden mean 
7, and the other one, to V^— 1. We also have on the crit- 
ical attractor, cycles (periodic and nonperiodic) related 
to the previous ones by symmetries, and therefore have 
identical properties. We can consider them as artefacts 
of the definition of the RG-transformation (for instance, 
when the frequency is equal to the golden mean, the RG- 
transformation can be modified such that the cycle of 
period three H,^ becomes a fixed point). 
The critical attractor is not irreducible [ P7|J2q| because it 
contains all the fixed points, periodic cycles etc., cor- 
responding to specific frequencies (e.g., quadratic fre- 
quencies). For a typical frequency, we expect the RG- 
trajectory to visit a subset which is dense in the whole 
attractor. 

Two invariant tori with frequencies lui and lu2 belong 
to the same universality class in the following sense : The 
properties of these tori at the threshold of their breakup 
are given by the analysis of a single hyperbolic attractor 
of a given renormalization transformation. These two tori 
display the same set of scaling factors and critical expo- 
nents (eigenvalues of the linearized map at each point of 
the attractor), with a different probability distribution, 
depending on the distribution of the entries in the con- 
tinued fraction of the frequency. 

Based on our numerical results, we can speak of a sin- 
gle generic universality class for the breakup of invariant 
tori for Hamiltonian systems with two degrees of free- 
dom. The universal properties are obtained from a single 
critical attractor of a RG-transformation. The critical at- 
tractor itself is not universal as it depends on the specific 
implementation of the transformation, but the properties 
derived from it, are universal. 



IV. FRACTAL DIAGRAM 

In order to visualize the critical surface S in the en- 
larged space, in which the critical attractor is contained, 
we represent a section of S by computing the fractal 
diagram Sc{tu) of the following one-parameter family of 
Hamiltonians : 

1 



7^"i?e ^ oo for £ > ec(7), 
WH^^Ho fore<ec(7). 



(4.2) 
(4.3) 



H,iA,ip)=u;o-A+^{n-Af 

+e (cos ifi + C0s((y3i + IP2)) 



(4.1) 



We compute ec(w), defined by Eqs. (4.2)-(4.3), with 
TZ the renormalization transformation constructed for 
quadratic Hamiltonians (p]q). We base our analysis on 
the conjecture that this £c{^) coincides with the thresh- 
old at which the torus with frequency to breaks up. Each 
frequency lo of the diagram is chosen at random in the 
interval [0.5,1]. In order to iterate TZ (defined by a se- 
quence of operators TZa^ ) , we compute the ten first entries 
in the continued fraction expansion of w, and the tail of 
this expansion is filled with a sequence of 1. The cut-off 
parameter for TZ (the restriction of the Fourier series to 
the modes v such that max^ \vi\ < L) is L ~ \Q . 
The curve £c{uj) depicted on Fig. |j, indicates the fractal 
structure of the critical surface S =_U„'5(w). For ratio- 
nal frequencies, £c is zero. In Fig. y, we observe strong 
resonances (Arnold tongues) close to rational frequencies 
^ = pI 1 with small q. Similar pictures have been ob- 
tained for the standard map [29|-|3^. 
The curve ec(w) for uj G [0, 0.5 [ is obtained from the one 
for ui g [0.5, 1] by the symmetry eduj) = £c(l ~ ^) which 
can be found by applying the canonical transformation 

{A,^)^{A\v') = {fA,T^), 

where T denotes the transposed matrix of 

l^ 1 

The effect of this transformation on H^ is just to change 
the frequency from lo to 1 — u. 

From Fig. 0, we can see a numerical indication that the 
golden mean torus (with frequency 7 « 0.618) is the 
most robust one for the one-parameter family ( |4.l| ) : 
£c{l) > ec(w) for all other uj G [0.5, 1]. 

For each line going through the origin in the MP- 
parameter plane for Hamiltonians 

H,{A,ip)^u:o-A + ^{n-Af 

-t-Af cos(^i + Pcos((^i + (P2), 

we can calculate in the same way, the critical coupling as 
a function of the frequency lo. This leads to pictures sim- 
ilar to Fig. ^. Putting together all these sections, we ob- 
tain a critical surface Sc{M^ P, to) for this two-parameter 
family of Hamiltonians and tori of frequencies to. The 
universal critical attractor is contained in a surface of 
this kind edifu },lo), but with an infinite number of 
AJ) 



ered space of Hamiltonians. 



where J7 = (1, 0) and a;o = {^, ^1), where lu is the fre- parameters /^r which are the coordinates of the consid- 

quency of the invariant torus. 

For the golden mean 7, we checked that £0(7) coincides 

with the value of the coupling where the RG-iteration 

starts to diverge, i.e., as n tends to infinity, we have 



V. A DIMENSIONAL RENORMALIZATION 
SCHEME 

In this section, we construct an approximate renor- 
malization scheme based on simple dimensional argu- 
ments. It is close to the type of transformations con- 
sidered by Escande and Doveil [p4|-p7[, and by MacKay 
and Stark [g6[. This scheme is built by arguments based 
on the dimensional analysis of the renormalization trans- 
formation. The aim is to see that the qualitative features 
of the renormalization flow obtained in the previous sec- 
tion, are already contained in a very simple scheme. 

We suppose that the initial Hamiltonian contains only 
the two Fourier modes Me**^^ and Pe*'^^ . The main scale 
is then determined by i/-i = (0, 1) and v>o = (1, 0). The 
next smaller scale is represented by the mode Vq = (1; 0) 
together with the next resonance Ui ~ (ao,l), where 
ao denotes the integer part of lo^^. The transforma- 
tion is a change of coordinates that eliminates the mode 
iz-i ~ (0,1), and produce the mode vi = (ao,l). As 
Ui = aQi/Q + u^i, the amplitude of this mode is MP'^° 
(to the lowest order). Then we shift the Fourier modes : 
the mode Vi (resp. Vq) becomes the mode i^o (resp. i^-i). 
Consequently the frequency uj is changed into co' accord- 
ing to the Gauss map (2.3). 



We obtain thus the following RG-scheme : 

M' = kiP, 
P' ^k2MP"-°, 

Lo' = [uj^^] - ao. 



(5.1) 
(5.2) 
(5.3) 



where Oq = [w^^]. In general, ki are f unct io ns o f oj, so the 
RG-scheme is equivalent to a system (5T)-(5^) driven by 
the Gauss map. 

To give an example, we consider ki — k2 — Lu^'^{ao + 
a)^ , where a is determined by the inverse of the Gauss 
map a' = l/(ao -I- a). We have chosen these coefficients 
identical to the rescaling coefficient of the constant term 
in the actions for the renormalization explained in Sec. O. 
This choice is to a large extent arbitrary, and is meant 
only as an illustration. The idea is that any function 
that reflects the Gauss map should lead to qualitatively 
similar results for the attractor and to the same critical 
exponents. 
We denote TZa,, the following map, where ao = [(^^^] ■ 



7^ao ■.{M,P,Lu,a) 



{M',P\uj',a'). 



For a given frequency lu — [ai,a2, . . .], the transforma- 
tion TZ is a sequence of TZai ■ This transformation TZ 
has two main domains : one where the iteration con- 
verges to Af = P = (trivial fixed set), and one where 
the iteration diverges. If the frequency is quadratic 
w = [&i,fe2, • • ■ ,6t, (ai,a2, . . . ,as)oo], there is a hyper- 
bolic (nontrivial) fixed point for the transformation TZg — 
TZa,TZas-i • • • ^oi • The critical surface is the codimension 
1 stable manifold of this nontrivial fixed point. 



For nonquadratic irrational frequencies, the critical sur- 
face is the codimension 1 stable manifold of a strange 
chaotic attractor : the chaoticity comes from the random 
sequence of the TZa^- We depict this critical attractor 
on Fig. o, computed using a frequency whose continued 
fraction expansion is a random sequence of 1 and 2 with 
p = 1/2. The largest Lyapunov exponent pS-WG] k gives 
the link between various Hamiltonians (of a given one- 
parameter family) near the critical surface. It measures 
how far are two Hamiltonians Hi et H2 (near the critical 
surface) , as we iterate renormalization TZ : 



7^"i/l - TZ"'H2 



\Hi - H2). 



We numerically find that the value given by this simple 
scheme, is in good agreement with the one given by the 
complete RG-transformation, k « 0.68 for p = 1/2. This 
exponent depends on p because it is computed for a given 
RG-trajectory which visits the different regions of the 
attractor with some p-dependent distribution. For p = 1 
(golden mean), k ~ 0.49 and for p = 0, k ~ 0.89. 
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APPENDIX: EQUATIONS DEFINING THE KAM 
TRANSFORMATIONS FOR 



HAMILTONIANS (2.9) 



In this appendix, we describe one step of the eli mina - 
tion procedure iJ k-> iJ', defined for Hamiltonians (p!s|), 
by considering that f^^' depends on a small parameter e, 
such that I^/(^^ is of order 0(e), where I^ f^J'> denotes 
the nonresonant part of f^-'\ i.e.. 



r/^^'(<P)= Y] f^^e^""" 



E 

vei- 



where fj^-'' denotes the Fourier coefficient of Z*^^-* with 
frequency vector 1/. We define Ho as 



iIo(A)=u;o-A + (/(2))(n.A)2 



(Al) 



In order to eliminate the nonresonant modes of /'^^■' to 
the first order in £, we perform a Lie transformation U : 
{A, if) i-^ (A', if') generated by a function S of the form 



S{A,ip) = i^Y'^^\ip){n ■ Ay + an ■ if. (A2) 



^E 



The expression of th e Hamiltonian in the new coordinates 
is given by Eq. (2.7). The first order in the perturbation 



of this Hamiltonian is V + {S, Hq}. Then S is determined 
by the foUowing condition 



-{S,Ho} + I-V^O. 



(A3) 



The constant a ehniinates the hnear term in the (J7 • A)- 
variable, (/(i)), by requiring that {{S,Ho}) + {f'^^^)n ■ 
A = 0: 



(/(^^ 



2f]2(/(2))' 

and Y^^") is determined by 
ia;o-5y(°)+r/("^ ^ const, 



0, 



(A4) 



(A5) 
(A6) 



for j > 1 , where d denotes the derivative with respect to 
the angles : d = d/dip. These equations are solved by 
representing them in the Fourier space : 



y''\^)= E 



/, 



(0) 



^IVip 



uei- 



Y^^\ip)= E 



vei- 



ujn ■ V 



u;o v 



/(^■)-2(/(2))fi.i.rj^--i: 



^iv-^p 



for j > 1. Then we compute H' = H oU he calculating 



recursively the Poisson brackets S''H — SS'' ^ 
k > 1. Denoting H^ = S'^H, H' becomes 



^' = E 



fc=0 



IT' 



We expand H' in power series in the actions 

oc 

H'{A, vp) = u;o • A + E f'^'^ i^^ ' ^y 

3=0 



H, for 



(A7) 



(A8) 



The Hamiltonian H' is expressed by the image of the 
functions f^^^ given by the following expressions 



/ 



'(j) 



oo fU) 



fc=0 



fc! 



(A9) 



where 



necessary to compute it for j too large as its contribution 
in H' might exceed the truncation in the actions. More 

precisely, we compute /^ for j = 0, . . . , J + fc( J — 1) 
if k is lower than fco = [K/J] + 1 (where [x] denotes the 
integer part of x), and for j — 0,...,J + K— kiik > ko- 
For instance, if we truncate at J = 3, we compute SH 
up to order (Jl • A)^. 



1=0 



U) 



r^ 



1=0 



for fc > 1 and j > 0. Numerically, we co mpu te /'^■'■' for 
j = 0, 1, . . . , J, by truncating the series (A9) to a finite 
sum over < k < K. For the calculation of fi^ , it is not 
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FIG. 1. Projection on the plane {f^Wfi]^), of the critical FIG. 3. Values of the rescalinss {Xi,Xi+i) after i iterations 

attractor for a frequency whose continued fraction expansion '^^ the critical attractor ot h ig. |1|. 
is a random sequence of 1 and 2 with probability P(l) — 1/2. 
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FIG. 2. Projection on the plane {fui , f^i), of the critical 
attractor for a frequency whose continued fraction expansion 
is a random sequence of 1 and 2 with probability P(l) = 9/10. 



FIG. 4. Fractal diagram edoj) for the one-parameter family 
of Hamiltonians ( [l.l| ) . 
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FIG. 5. Critical attractor for the approximate renormal- 
ization scheme for a frequency whose continued fraction ex- 
pansion is a random sequence of 1 and 2 with probability 
P(l) = 1/2. 
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